Scaling picture of magnetism formation in the anomalous f-systems: interplay of the 

Kondo effect and spin dynamics 
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Abstract 



Formation of magnetically ordered state in the Kondo 
lattices is treated within the degenerate s — / exchange 
and Coqblin-Schrieffer models. The Kondo renormal- 
izations of the effective coupling parameter, magnetic 
moment and spin excitation frequencies are calculated 
within perturbation theory. The results of one-loop scal- 
ing consideration of the magnetic state in Kondo lat- 
tices are analyzed. The dependence of the critical val- 
ues of the bare model parameters on the type of the 
magnetic phase and space dimensionality is investigated. 
Rcnormalization of the effective Kondo temperature by 
the interatomic exchange interactions is calculated. An 
important role of the character of spin dynamics (ex- 
istence of well-defined magnon excitations, presence of 
magnetic anisotropy etc.) is demonstrated. The regime 
of strongly suppressed magnetic moments, which corre- 
sponds to magnetic heavy-fermion system, may occur 
in a rather narrow parameter region only. At the same 
time, in the magnetically ordered phases the renormal- 
ized Kondo temperature depends weakly on the bare 
coupling parameter in some interval. The critical be- 
havior, corresponding to the magnetic transition with 
changing the bare s — f coupling parameter, is inves- 
tigated. In the vicinity of the strong coupling regime, 
the spectrum of the Bose excitations becomes softened. 
Thus on the borderline of magnetic instability the Fermi- 
liquid picture is violated in some temerature interval due 
to scattering of electrons by these bosons. This may 
explain the fact that a non-Fermi-liquid behavior often 
takes place in the heavy-fermion systems near the onset 
of magnetic ordering. 

75.30.Mb, 71.28-fd 



1. Introduction 

Anomalous 4/- and 5/-compounds, including so-called 
Kondo lattices and heavy-fermion systems, aire, studied 
extensively starting from the middle of SOaJu. From 
very beginning of these investigations it became clear 
that the effects, connected with regular arrangement of 
the Kondo centres (rare-earth or actinide ions) play a 
crucial role in the physics of such systems. When passing 
from one magnetic centre to the Kondo lattice, two main 



new features appear. First, provided that the strong cou- 
pling regime takes place, the Abrikosov-Suhl resonance 
in the one-site t-matrix leads to formation of a compli- 
cated band structure near Ep on the new energy scale 
(the Kondo temperature Tk) Wfr^ sharp peaks and pseu- 
dogaps in the density of state^u. This provides a com- 
mon explanation of the heavy-fermion behavior. Sec- 
ond, the competition between the Kondo screening of 
magnetic moments and-.ij||J.ersite magnetic interactions 
has a great imnortanceErElu^Ll. Following to the old pa- 
per by DoniachB, it was believed in early works that this 
competition leads to total suppression of either magnetic 
moments or the Kondo anomalies. However, more recent 
experimental data and careful theoretical investigations 
made clear that the Kondo lattices as a rule demonstrate 
magnetic ordering or are close to this. This concept was 
consistjfintly formulated and justified in a series of our 
papersLTtJ. A very important circumstance is that in- 
terspin coupling between the Kondo sites results in a 
smearing of singularities in electron and magnetic prop- 
erties on the scale of the characteristic spin-dynamics 
frequency lJ. At the same time, ZU itself acquires renor- 
malizations resulting in its decrease due to the Kondo 
screening. A simple scaling consideration of this rcnor- 
malization process in the s — f exchange modelEj yields, 
depending on the values of bare parameters, both the 
"usual" states (a non-magnetic Kondo lattice or a mag- 
net with weak Kondo contributions) and the peculiar 
magnetic Kondo-lattice state. In the latter state, small 
variations of parameters result in strong changes of the 
ground-state moment. Thus a characteristic feature of 
heavy fermion magnets — high sensitivity of the ground- 
state moment to external factors like pressure and dop- 
ing by a small amount of impurities — is naturally ex- 
plained. At the same time, only the simplest s — f model 
was considered in Ref£3, and the equations obtained 
were not investigated in detail. Therefore a number of 
important features of the Kondo magnets were not de- 
scribed. 

Recently, a num- 

ber of anomalous /-systems (UxYi_a;Pd3, VPt^-xPdx, 
UCus-xPd^:, CeCug-ajAux, UxThi_a;Bei3 etc.) demon- 
strating the so-called non-Fermi-liquid (NFL) behavior 
have become a subject of great interest (see, e.g., the 
reviewsliJ). It should be noted that such a behavior is 
observed not only in alloys, but-jalso in some stoichiomet- 
ric compounds, e.g., CerNiaEl, CeCu2Si2,CeNi2Ge2y. 
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These systems possess unusual logarithmic or power-law 
temperature dependences of electron and magnetic prop- 
erties. It is a common practice to discuss such a behavior 
within the one-impurity two-channel Kondo modelO'ES. 
However, the NFL behavior is typical for systems ly- 
ing on the boundary of magnetic-, ordering and demon- 
strating strong spin fluctuationaiJ. So, many-center ef- 
fects should play an important role in this phenomenon. 
At the same time, for a number anomalous /-systems 
as Sm3Se4, Yb4As3, as well as for the only "moder- 
ately heavy-fermion" d-system Yi_a;Sc2:Mn2, the heavy- 
fermion state itself seems to be connected with peculiari- 
ties of intersite couplings (e.g., frustcatintis), rather than 
with the one-impurity Kondo effectO'UO. Thus the in- 
terplay of the Kondo effect and intersite spin dynamics 
results in very rich and complicated picture rather than 
in trivial mutual suppression. 

The aim of the present paper is a systematic investi- 
gation of formation of the magnetic Kondo-lattice state 
and of its properties for various magnetic phases depend- 
ing on the character of spin dynamics. 

In Sect. 2 we introduce main theoretical models which 
enable one to treat the Kondo effect in a lattice in differ- 
ent cases with account of orbital degeneracy. In Sect. 3 
we analyze the properties of the localized spin subsys- 
tem for these models in the absence of s — / interac- 
tion. The Kondo renormalizations in the paramagnetic 
state with account of spin dynamics are considered in 
Sect. 4. In Sects. 5 and 6 we calculate the Kondo cor- 
rections to the electron spectrum (and thereby to the 
effective s — f coupling) and magnon frequencies in the 
magnetically ordered phases. In Sect. 7 we write down 
the scaling equations for the effective s — f parameter 
and spin-fluctuation frequency which plays the role of 
a cutoff for the Kondo divergences in concentrated /- 
systems. In Sect. 8 we generalize the scaling approach 
for the ordered state by including the renormalization of 
the residue of the spin Green's function at the magnon 
pole. The most simple large-iV limit in the Coqblin- 
Schrieffer model where spin dynamics is unrenormalized 
is considered in Sect. 9. The scaling picture for finite N 
is discussed in Sect. 10. In Sect. 11 we discuss the critical 
behavior near the magnetic phase transition and discuss 
a possibility of the Fermi-liquid picture violation. The 
scaling behavior in the s — f exchange model with a 
large orbital degeneracy is investigated in Sect. 12, and 
an explicit description of the non-Fermi-liquid behavior 
is obtained in this limiting case. 



2. Theoretical models 

To treat the Kondo effect in a lattice we use the s — d{f) 
exchange Hamiltonian 

^ = E^kcLck. + i?/+i?s/ = i/o + ifs/ (1) 

where tk is the band energy. We consider the pure spin 
s — d{f) exchange model with 

i?/ = E JqS-qSq, Hsf=- ^ 4k' Sk-k'C^Q/jcj^QCk'/S 
q kk'a^ 

(2) 

where Si and Sq are spin operators and their Fourier 
transforms, a are the Pauli matrices. For the sake of 
convenient constructing perturbation theory, we explic- 
itly include the Heisenberg exchange interaction with the 
parameters Jq in the Hamiltonian, although in fact this 
interaction can be the indirect RKKY coupling. Ex- 
panding the s — d{f) interaction in spherical functions 
we have 

/kk' = Yl ^'^/-(^k, 0k)i1™(^k', <^k') (3) 

Im 

Hereafter we retain in (^) only one term Ii = I {I — 2 
for c?-electrons and ^ = 3 for /-electrons). Introducing 
the operators 

Ckm<T = i'(47r)i/2^ka>";m(6'k, 0k) (4) 

which satisfy, after averaging over the angles of the vec- 
tor k, the Fermi commutation relations, we reduce Hgf 
to the form 

kk'mQ/3 

Assuming the electron and spin excitation spectrum to 
be isotropic, in final expressions for self-energies we can 
perform averaging over the angles of all the wavevectors 
and use the orthogonality relation 

sin9dedHYilie,(b)Yi^'i0,(j)) = Smm' 

Thus the factors of [I] = (2Z -I- 1) occur in any order of 
perturbation theory, and we have to replace /" [Z]/" 
in the "connected" terms of perturbation expansion in 
comparison with the "standard" s — d model {1 = 0). 

It is worthwhile to remembtpjpain results for the one- 
impurity version of this modelHj'Ell. Perturbation theory 
treatment leads to occurrence of infrared divergences. 
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Provided that / < 0, the characteristic energy scale (the 
Kondo temperature) 

Tk = De^j>{l/2Ip) (6) 

occurs, where D is of order of bandwidth, p is the 
bare density of electron states at the Fermi level with 
one spin projection. At T ~ Tk the effective s — / 
interaction becomes very large and the system enters 
the strong coupling regime. The electron energy spec- 
trum in this region is determined by the presence of the 
Abrikosov-Suhl resonance of the width Tk- Properties of 
the ground state and character of the low-temperature 
behavior depend crucially on the parameters S and [I]. 
At 25* = \l] the Fermi-liquid singlet state occurs. At 
2S > [I] the localized moment and logarithmic behavior 
of electronic characterisitics retains, but the replacement 
S ^ S — [l]/2 takes place. At 25* < [I] a very interesting 
"overcompensated" regime occurs. Last time, the par- 
ticular case of this regime with 5 = 1/2, [/] = 2 (the two- 
channel Kondo m|©del describing the non-Fermi-liquid 
behavior, see, e.g.JlZllla) is a subject of great interest. 

In the case of the periodic model presence or absence of 
the strong coupling regime depends also on the character 
of intcrsitc spin-spin interactions which are described by 
Hf. This factor will be analyzed in detail below. 

The s — d{f) model does not takejutij-apiount scat- 
tering by orbital degrees of freedomEScJ^EJ^cJ. Another 
important model, which is used frequently to describe 
the Kondo effect, is the Coqblin-Schrieffer model. For 
its periodic version with the f — f exchange interaction 
of the SU{N) form we have 

Hf-\T.J^ E X^'X^'^, (7) 

q M,M' = -S 

U ^ T \^ vMM'J 

where AT*^*^ are the Fourier transforms of the Hub- 
bard's operators for the localized spin system 

XMM' ^ \iM){iM'\ 

S = {N — l)/2 is the total angular momentum (this 
notation is used for the sake of convenience and has a 
somewhat different meaning in comparison with the spin 
S in the s — f model) , the operators 

CkM = ^'iallm^^-'no- (8) 

ni(7 

possess, after averaging over the angles, the Fermi prop- 
erties due to the orthogonality relations for the Clebsh- 
Gordan coefficients C. As well as for the model (0), 



we will assume that this averaging should be performed 
when calculating the Green's functions. The Hamilto- 
nian Hsf with J < can be derived from the degenerate 
Anderson-lattice model for rare-earth compounds (LS- 
coupling) ; a Hamiltonian of the same farm occurs in the 
case of jj-coupling (actinide systems )E3. 

For = 1/2,7V = 2, Z = the models (|) and (0) 
reduce to the standard s — / model with 5 = 1/2 and 
coincide. The ground state in the one-impurity Goqblin- 
Schrieffer model is similar to that in the s ~ f model 
with 25 — [I], i.e. complete screening of the localized 
moment and a Fermi-liquid picture take place. Due to 
another structure of perturbation theory for the Goqblin- 
Schrieffer model, we have to replace 2 ^ in (^. Thus 
the role of the degeneracy factors in both the models 
under consideration is different: the expression for Tk 
does not contain the factor of [I] = (21-1-1) in the model 
(H), but contains the factor of N in the model (0). Pecu- 
liarities of the Coqblin-Schrieffer model are determined 
by that the transitions between any values of localized 
/-state projection M are possible, so that the number of 
excitations branches is large. We shall see that this may 
result in essential modifications of magnetic behavior. 

The interaction Hf in (|^) can be obtained as an indi- 
rect RKKY-type interaction which arises to second order 
in Hsf. Using the standard Heisenberg interaction (as 
in (|2|)), where only M —f M ± I transitions are allowed, 
is incovenient since other transitions acquire an energy 
gap. However, inclusion of this interaction does not lead 
to a strong change of the physical picture. The standard 
angular momentum operators on a site are expressed in 
terms of the A-operators as 

S+ = ^(5 - Af)i/2(5 + M + i)i/2xf +1'^^ (9) 

M 
M 

Of course, both the s — f exchange model and Coqblin- 
Schrieffer model are some idealizations of realistic situ- 
ation. Choice of adequate model for a given compound 
depends mainly on the relation between the width of 
/-level and the spin-orbital coupling parameter. If the 
broadening of /-level due to either the hybridization or 
direct / — / overlap is larger than its spin-orbital split- 
ting we should consider the latter after the transition 
from atom-like /-states to the crystal states. Usualfet 
the orbital moment is quenched in Bloch-like statesE£l 
and therefore only the spin moment should be taken into 
account when considering the interaction with conduc- 
tion electrons. The Coqblin-Schrieffer model was ini- 
tially proposed to describe cerium and ytterbium sys- 
tems, especially diluted onesEHl. However, as it is clear 
now, the situation for cerium compounds is more compli- 
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cated. First, the spin-orbit coupling for cerium is in fact 
not too large (about 0.25 eV) and comparable with the 
width of the virtual /-level. Second, a number of cerium 
system, including the pure a— cerium have rather large 
/ — / overlap (the relative role of / — / overlap-iand hy- 
bridization is discussed in detail in the reviewcZI). Thus 
the applicability of one of these models should be con- 
sidered separately for any specific compound. We shall 
demonstrate below that results of the scaling considera- 
tion for the Coqblin-Schrieffer model and s — / exchange 
model are essentially different. 

The model (|) may be generalized to include two mag- 
netic /-configurations with the angular moments J and 
J' so that 



H. 



J' fi,jm^ J' fi,jm 



vMM' J 
'^k'-kCk',„'Ckn 



(10) 



(we restrict ourselves for simplicity to the case of 
jj-coupling bearing in mind uranium compounds). 
However J— tijii Hamiltonian has a complicated tensor 
structura23'Ea and does not enable one to calculate an 
unique energy scale by using perturbation expansion. 
Such an energy scale can be obtained starting from the 
low-temperature regime and read£3 



Tk = D exp 



-1 Ip 



2J+1 
2 J' + 1 



1 



(11) 



(note that the exponent in ( [III ) for the case J' — differs 
by an unity from the correct result; such a difference is 
typical for the mc.thpds which are in fact based on the 
large- expansionE3c3) . In the case J > J' the situation 
for the model (^l]) is similar to that for the s — / model 

with 25* > m. 



3. Ground state and spectrum of spin excitations 

In the ground ferromagnetic (FM) state to zeroth ap- 
proximation in / we have (S'f ) = S, and the spin-wave 
spectrum for the model (0) reads 



,FM 



(5) = 25( Jq - Jo) 



(12) 



For the model (|7|) , we also assume the magnon character 
of spin dynamics. Then the spin-wave spectrum can be 
found by linearizing the equation of motion for the "spin 
deviation" operator ^:i^q^. There exist 25* = iV — 1 spin- 
wave modes with the frequency Uq — ajq*^(l/2), which 
correspond to the transitions S* — > M < 5. 

Now we discuss an antiferromagnet which has the spi- 
ral structure along the a;-axis with the wavevector Q 

( Jq — Jmin) 

{S^^) = ^cosQR,, {S^) = ^sinQR,, = 



We introduce the local coordinate system 

S-f = 5f cos QR, - Sf sinQR,,, 

Sf = Sf cos QR, + S-f sin QR,, Sf = Sf 



(13) 



Hereafter we consider for simplicity a two-sublattice 
AFM (2Q is equal to a reciprocal lattice vector, so that 
cos^ QRj = 1, sin^ QRj — 0). Passing to the spin devi- 
ation operator we derive in the spin-wave region 

Hf = const + Y^iCqb^bq + ^Dqib^qbq + b^bl^)] (14) 
q 

Cq = S'(jQ+q + Jq — 2Jq), I?q = S'(Jq — Jq+q) 

Diagonalizing (llj) we obtain the spin-wave spectrum 



.AFM 



iS) = K-Diy/' = 25(Jq-jQ)V2(jQ^^_ Jq)1 

(15) 



/2 



The corresponding transformation to the local coordi- 
nate system for the model (^ has the form 



Then we obtain 



^*^'(l-fcos QRJ+Xr^^-*^'(l-cos QR,)] 



q M.M' = ~S 

(16) 



^^/ = -^E E [(^. 

kq M,M'=-S 



M,M' I Y-M-M'\ 
q + -^q ) 



M,M' 

q+Q 



X 



-M,-M' 

q+Q 



)]4M'Ck- 



qAf 



(17) 



with 



7(1,2) 

^q 



In the mean field approximation we have 
S 

M=-S 



(18) 



with um = (X^^). The usual AFM state turns out 
to be unstable for simple lattices provided that only 
the nearest-neighbor interaction is taken into account 

(2) 

(Jq =0). Indeed, in this case any state with 



Y "Af 1, X! ^ ^ 



M>Q 



M<0 
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has the same energy {Hf) = 0. When introducing the 

(2) 

next-neighbor ferromagnetic interaction (Jg < 0) the 
AFM state with 

ns = 1, UAKS = (19) 

is stabiHzed. This state will be referred to as AFMl. 
The corresponding "spin-wave" spectrum contains N 
branches. The mode, which corresponds to the transi- 
tion S^-S, has the frequency = = 1/2). 
Other N — 2 modes have a ferromagnetic type and pos- 
sess the energy 

1/2, Jq 



7(2) 



'q ^ ~q ''q 
Provided that next-neighbor interaction is antiferro- 

(2) 

magnetic too ( Jq > 0, this case is referred to as AFM2), 
minimization of ( |l8| ) yields 

2/N = 2/(25+1), N even 

Jo/iSJo + 4^'^), N odd 
(20) 



riMKO = 0, nM>o = K ^ 



For odd N we have 



(2) 



(2)n 



no = Jr/('5Jo + JD - Jo/iSJo) « if (21) 

Then, according to (^) , for iV > 2 the sublattice magne- 
tization turns out to be reduced in comparison with S 
already in the mean-field approximation: 



5 + 1/2,5 half-integer 
(5+ 1)(1 - no), 5 integer 



(22) 



Then there exist N'^/A {N even) or {N - 1)^/4 {N 
odd) "antiferromagnetic" modes which correspond to the 
transitions from M > to M' < and have the fre- 



(a) 

quency ujq 



,AFM ( 



Besides that. 



for odd 
o 



q {S ^K) 

N there exist peculiar modes with the frequencies tjq, 
which describe the transitions M and are deter- 
mined by the equations 

(23) 

(the plus sign corresponds to Af > and the minus sign 
to M < 0), so that to lowest order in no 



Tno {€% - Dl)/C^ 



(24) 



^ql ^-"■^qj'^ql 

Thus one of the solutions to (Esl) describes the mode 

which is very soft for Jq <C 5 Jo. Since the Kondo sin- 
gular contributions (see below), which are cut at this 
mode, are large, this may result in a tendency to the 
destruction of magnetism. 

Really, the considered ground states and excitalion 
spectra are strongly influenced by the crystal fieldn3't3'Lil. 
In particular, large gaps and degeneracy lift may occur in 
the above-discussed "additional" modes. This question 
needs a special consideration. 



4. Renormalization in the paramagnetic phase 

The Kondo-lattice problem in the paramagnetic state 
describes the process of screening of localized magnetic 
moments. 

To find the renormalization of the effective s — f ex- 
change parameter we consider the electron self-energy. 
In the second order in / we have 



sf(£;) = /2pV 



1 



^ E- 
q 



where 



^^1^5(5 + 1), model (2) 
\ l~l/N^, model (0) 



(25) 



(26) 



To construct a self-consistent theory of Kondo lattices 
we have to calculate the third-order Kondo correction 
to the self energy with account of spin dynamics. Such 
calculations were performed in RcfsJj within the simplest 
s — / model. In our case we obtain for the singular 
contribution (which contains logarithmic divergence in 
the absence of spin dynamics) 



dc^y^ Jq(cj) 

q,p 



(27) 



where TV = 2 in the model (g), J7q(ti^) is the spectral 
density of the spin Green's function for the Hamiltonian 
Hf, which is normalized to unity, the spin dynamics is 
neglected in the denominator, which is not connected 
with the Fermi function nk = n(tk)- The Kondo renor- 
malization of the s — / parameter / I^f ~ I + Sl^f 

f3l (2) 

is determined by "including" IniEj^ {E) into ImSj^. {E) 
and is given by 



= 



SI. 



(28) 



(here and hereafter we have to put E = Ep = 0, k = kp 
in the expressions for lef)- To concretize the form of 
spin dynamics in the paramagnetic phase we use below 
the spin diffusion approximation 



(29) 



{V is the spin diffusion constant) which is correct at 
small q,^; and is reasonable in a general case. Then we 
derive 



"J 2 E- ik-q - iVq^ 

q ^ 



(30) 



1 



E-tu- 
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To calculate the correction to the effective magnetic 
moment we treat the static magnetic susceptibility 

X = {S\S') = [ ' dA(exp(AiI)5'^exp(-AiJ)S'^) 

(31) 

Expanding to second order in / we derive (cf. Ref.E 



(32) 



where we have introduced the notation 

R = 



_ ( [I], model (2) 
\N/2, model (0) 



(33) 



The spin-fluctuation frequency in the paramagnetic 
phase is determined from the second moment of the spin 
Green's function 



(S'-q, Sq)/ ('5'^q, <5'q) 



To second order in / we derive (cfBi ) 
(c.2)o = l5(5 + l)^(Jq_p- Jp)2 



(34) 



(35) 



(g). For a ferromagnet the electron spectrum possesses 
the spin splitting, Ei^„ = ik ~ <7[l]IS. The second-order 
correction to lef is determined by the corresponding elec- 
tron self-energies: 



SI 



j:[IHe)]/{2S[i]) 



(39) 



with 



E — tk-q 



,FM ■ 



^^l'iE)^2RI^Sj2 



E — tk-c 



,FM 



(40) 



For an antiferromagnet the electron spectrum contains 
the AFM gap, 



Ek — -(tk + ik+q) ± [^(ik - ik+q) 



{[l]ISf]^/^ 
(41) 



The renormalization of / is obtained from the second- 
order correction to the anomalous Green's function 



[l]IS- 



Si^.^yq(i^) 



so that 



(^^q/Wq = (1 - "q)'5 '5'^// S ef = "(1 " aq)L 



(36) 



where he have taken into account spin dynamics by anal- 
ogy with (HI). Passing into real space yields for the 
quantity 



R 



sin kpR 
kpR 



Slef = -S 



AFM 
k,k+q 



{E)/{S[l]) 



(42) 



The calculation of the off-diagonal self-energy gives 

nk-q(S - <k-q) 



{E - tk-q)2 - 



[1-cosqR]/^ J^[l-cosqR] 



(43) 



R 



(37) 



In the approximation of nearest neighbors at the distance 
d, a does not depend on q: 



sin kpd 



(38) 



so that we may use a single renormalization parameter, 
rather than the whole function of q. 



5. Effective interaction in magnetically ordered 
phases 

Now we investigate the renormalization of the s — f inter- 
action in FM and AFM phases. First we treat the model 



To calculate the corrections to lef in the Coqblin- 
Schrieffer model we consider the Green's function 



(((CkAf |Ck+qM>>S ~ ((Ck,-Af l4 



k+q,-M 



))e) 



M>0 



(44) 



which determines the "magnetization" (FM case, Q = 0) 
or "staggered magnetization" (AFM case) of conduction 
electrons. 

For a ferromagnet the mean-field electron spectrum 
reads -EkAf = ^k — I^ms- Calculating the second-order 
corrections we derive 

Slef = -2[{N m^l'iE) ^[t'iE)]/N (45) 

(remember that in the Coqblin-Schrieffer model the self- 
energies should be substituted at S' = 1/2). 
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In the AFMl state the mean-field electron spectrum 
for M = S,—S is given by (^) with the replacement 
tk t]s_ — I. For other M the spectrum is unrenor- 
malized. The renormalization of / in such a situation 
contains contributions of both FM and AFM types: 



5(J. 



tot 



jtot \ 



(f) AFM exchange contributions of the first order in 

'^q )+Sk,k+Q(£^> l/2Sl;^<icf^^d to the RKKY approximation 



(46) 



In the case AFM2 the electron spectrum is given by ( [4l| ) 
with S K/2, tk <k - IK. Besides that, for odd N 
there exists a branch of spectrum with M = 0, which is 
weakly renormalized due to smallness of | Jg^^ | . Then we 
obtain for even N 



nk - nk-c 



^ UJ + tl, 
k 



k-q 



(54) 



r T oy^AFM ( 77 AFM 



For odd N the contribution of the mode 



(47) 



the second summand in ( p4| ) being the w-dependent 
RKKY indirect exchange interaction. The function 
which determines the second-order corrections, is given 
by 



occurs: 



&IeS = -2 



^AFM 
^k,k+Q 



AFM 

q 



AFM 
pqw 



pqt^ 



2 nk(l - nu+p^q) + iV(±Wp)(nk 



"-k4 



k 



OJ - 



tk — ^k+p-q T '^p 



n+ 

k,k+Q 



where 



n-k-qiE — ik-q 



n^{E- 



ik-q) 



(49) 



6. Renormalizations of ordered moment and magnon 
frequency 

To investigate the magnon spectrum of an antiferromag- 
net in the model (||), we calculate the retarded Green's 
function of spin deviation operators in the local coordi- 
nate system 



i(^) = HKK))., rq(c.) 



''T'0^kk+QiE)]/K (48) where N{uj) is the Bose function (note that (/)+qj^ = 
— (/)pq_^), cjp is the magnon frequency to zeroth order 
in / and 1/25*, given by (|l5[), the terms in ( |53| ) which 
contain the spin-deviation correlation functions describe 
the magnon anharmonicity. The expression (|5^) is valid 
also for aJerromagnet (Q = 0) provided that q is not too 
small, cfS; such an approximation is sufficient to treat 
the Kondo divergences. 

We have to take into account singular s — f contribu- 
tions to the averages in (^) . These are due to the zero- 
point magnon damping and are obtained by using the 
spectral representation for the Green's functions (|5l|), 
( p2|) in the RKKY approximatioro. Since ImJq°* ^ w 
(\lo\ <C Ep) the corresponding integral over frequency 
contains logarithmic Kondo-like divergences smeared by 
spin dynamics (note that scattering processes correction 
to the damping described by the function (55) do not 
contribute the logarithmic terms). We derive 



(50) 



Writing down the equation of motion we derive (cf. the 
calculations for I = 0e3) 



rq(^) = 



(&-q&q> 



2.\niS>AFM , ^AFM 



'^^q+QOO^ 



q— u; 



{lo - C^^){lo + C^^^) + Dl 



qu; 



fq(^) 



qcj 



(t^-Cq^)(w + Cq_^)-|-L»2 



(51) 

(52) 



for an antiferromagnet and 



(57) 



(58) 



with 



where 



Q( jtot 



p 



[CpK 



FM 



k 



nk(l - nk-i-p-q) + iV(a;p)(nk - rik+p- 



(w -f ik - tkH 



p-q 



(59) 



AFM 
pOO 



pqu v^pqojj fQj. ferroVnagnet. The expressions 



8^ determine 



p 



[(2Jq + 2Jq_p - 2Jp - Jq+q - Jq)(6^6p) -al@gpt|SLg,iS^i]ilar correction to the (sublattice) magneti- 



zation 
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SS/S — ~—'S^S{b^b ) — — <i>''^*^''^^^''^ (60) According to (||), the magnetization of a ferromagnet 

S ^ *^ ^ is determined by 

Collecting all the singular s — d contributions to the pole ^^'^ ^ ^ Nn^, n_ — ^^("'^-q ) (^9) 
of (|5^) and taking into account the relation 

where the average in the right-hand side does not depend 

E7 ^ 7 ^ 72 / 7 \ , QyM iov M <S. Then we obtain 

"^p+k^^pqw — '^p+k-q'i^pOO ^ J \Jk-q/tk=tk=0 ZZ 

.gl^ ^ SS/S=-^Y.<o^^Su;^-Y^<'' (70) 

^ ' q 

which holds to logarithmic accuracy, we derive for the The sublattice magnetization in the AFMl case is 

singular correction given by 

Siu^fr ^ 2u;^Slo^^^' S/S=1-{N- 2)n_ - ^^^s (71) 

with 



q 

with R = [I]. In the case of a ferromagnet (Q =0) the ^ _ "S^ (X^^^ X^'^^) 

term (/)pq^ + (t)pqi^ (which is odd in lo) yields a contribu- ^-^ ^'^ 

tion to the pole of (|5^) and we have 

After cumbersome calculations we derive 

<S<*^(c.) = ~2RSY,{2Jp - 2Jq_p + Jq - Jo + ^/2^)a>^n^^^^^,^6_3) iV - 2 ^ ^fm^^fm ^ ^(/))_^ I'^o'^o*' ('^p 
p 2 

5<*^ = <5<*^(c.q) = -4i?5^(Jp + Jq-Jq_p- Jo)<i>^o^^ (64) (73) 



p 



This expression may be represented as N 

Scojcu^ = 2(1 - a^)6S/S (65) '5(c.(«))2 ^ A^(^^^'M)2(^afm ^ ^(a)) _^ 2i^C;q5<*^(c. 



with ^ith = in (|63D, (|62 

In the AFM2 case we have to put in (^ M < 0, M' > 



Y^Jr 1(6'"^)?,=^;^ [1 -cos qR]/^./R[l -COS qft] Then we have 

R R d ro\ f (5* + l/2)27i_, iV even 



(66) ^-(^)o [s{S+l)n^ + {S+l)Sno,N odd 

(75) 

For an antiferromagnet in the nearest-neighbor ap- 
proximation (Jq+q = — Jq) wc obtain from (62) where {S)o is given by (p2[), the average 



(JtJq/Wq = 55-/5 (67) - Z^(^-q ^q ) 

q 

Note that the results (|5l),(|6^) differ from those of RefJli does not depend on M, M' for M' > 0, M < 0, and Shq 

since only corrections arising from the static correlation the fluctuation correction to no- Restricting ourselves 

functions were taken into account in that paper. for simplicity to the case of even N, which corresponds 

The calculations of the magnon spectrum in the to a realistic situation for /-ions, we obtain 
Coqblin-Schrieffer model for FM and AFMl cases are 

performed in a similar way by calculating the Green's SS/S — — ^qop^^ (^p ^ ^p"'') (76) 

functions q 

Tq (W) = ((ATq |X_q ))c^,rq [uj) = {{X^ \X ^{a)) ±.^AF M , AF M ^a)^ 

(68) q ^ q ^ p p ^ ^ ' 

fThus the corrections to the sublattice magnetization and 

) , (pi) by the magnon frequency do not contain the factor of N in this 

in ^5^). (jase 
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7. Scaling equations 

Using the results of previous sections we can write down 
the system of scahng equations in the case of the Kondo 
lattice for various magaetic phases. We use the "poor 
man scaling" approachu. In this method one considers 
the dependence of effective (renormalized) model param- 
eters on the cutoff parameter C which occurs at picking 
out the Kondo singular terms. To find the equation for 

1) the con- 



For d = 2 the function rj^^^^ {x) vanishes discontinuosly 
at a; > 1, but a smooth contribution occurs for more 
realistic models of magnon spectrum. 

The singularities are retained if we introduce the gap 
in the spin-wave spectrum. For a ferromagnet is re- 



placed by ujo - 



Then we have, besides the singularity 



at X 



1, the second singularity at a; > 1: 



1, (40), 



lef we calculate in the sums in 
tribution of intermediate electron states near the Fermi 
level with C < tk+q < C + SC. Then we obtain for the 
model (|), and for the PM, FM and AFM2 phases in the 
Coqblin-Schrieffer model 



± 



LOq ■ 



2xUJe 



■In 



1 ±a; 



1 ± LUox/{uJo + LUex) 



(83) 



where ujex = maxwq. In the AFM case is replaced by 
uiq + and we derive 



Shf{C)^NplM~-^,)SC/C 



(78) 



~AFM 



(x) 



AFM 



^g(l 



[^2(i_^2)+^2Jl/2 



where ZU is a characteristic spin fluctuation energy, N = 2 
for the / model, r](x) is the scaling function which sat- 
isfies the condition ?7(0) = 1 Ithis guarantees the correct 
one-impurity limit, cf. Ref.EJ). For the paramagnetic 
phase we obtain 



(84) 



We shall see that the latter singularity may influ- 
ence considerably the behavior of the effective coupling 
constant. In the limit of strong magnetic anisotropy 
^o/i-^ex oo the singularity at a; —> 1 is very strong: 



For FM and AFM phases we have 



1 



(79) 



T]%x) = [i-x^y 



(85) 



r^'i7^)^T?iiN~l)m{j,)+mi7^)] (80) 



N 

U 

c 



^AFM^-^ = ((1 - U.l_^,/C^)-')t,=t,,=E. (81) 



C' 

Using the spin diffusion approximation ( p^ ) in ( [79[ ) and 
the approximations Wq^^ ^ ^ lu^-^^^ ^ q (which are 
justified, e.g., at small kp), we get 

x^"'" arctana; (82) 
±a;"Mn|l±a;| 

-x-'^\n\l~ x'^l, d = 3 
(l-x2)-i/26l(l-x2), d = 2 

where uJ = AVkp for a paramagnet, tU = W2feF for FM 
and AFM cases, 9{x) is the step function. One can see 
that the scaling functions for the ordered phases contain 
singularities at a; = 1. Presence of this singularity is a 
general property which does not depend on the spectrum 
model. The functions -qf^^x) and r)^^^\x) {d = 3) 

change their sign at a; = 2 and x — \/2 respectively. 





(x) 




(x) 




AFM 


(x) 



iO)/C 

When considering characteristics of localized-spin sub- 
system, the lowest-order Kondo corrections originate 
from double integrals over both electron and hole states 
(see (^H), (p9|), (p5|)). Then we have to introduce two 
cutoff parameters Ce and C/,. with Ce + Ch = C {C is 
the cutoff parameter for the electron- hole excitations). 
In the FM case for the Coqblin-Schrieffer model we have 
SCfi = —{N— l)5Ce due to the requirement of the 
number-of-particle conservation in electron-hole transi- 
tions (there exists iV — 1 "channels" for electrons and 
one "channel" for holes). For PM and AFM2 cases in 
the Coqblin-Schrieffer model, as well as for all the cases 
in the s — f model, the electron-hole symmetry is not 
violated and we have 50^ = —SCh- 

Taking into account (|3^, ([70|), ( [76| ) we obtain 

SSefiC)/S = VpSIefiC) = VNp^I^7j{-^)6C/C 

(86) 

where V = [I] for the s — f model, V" = 1 for PM and 
FM phases in the Coqblin-Schrieffer model and V = 2/N 
for the AFM2 phase. The renormalizations of spin-wave 
frequencies are obtained in a similar way from (|3p) , (|6S" 
(H), (Foh, (IttI), and are given by 



6LJ^f{C)/LJ = aSSef{C)/S = aVNp''rr]{--)SC/C 

(87) 
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where, in the nearest-neighbor approximation, 

( 1 - a PM 
a = < 2(1 - a) FM 
[ 1 AFM 

Introducing the dimcnsionlcss coupling constants 

5e/(C) = -NpI,f{C), g = -Npl 



((5+|5Zq)). = ^ + ((5+|5lq))r°' (93) 

(88) 

where the residue at the pole is determined by (K33): 



and replacing in the right-hand parts of (|78|), (|87]), ( |86| ) 

3 ^ gef{C), uj — > uJe/(C) we obtain the system of scaling 
equation 



dgef{C)/dC =-K 
dlnujefiC)/dC ^ aVA/N 
d\nSef{C)/dC ^VA/N 



(89) 
(90) 
(91) 



with 



A = A(C,cJe/(C)) = [.9e'/(C)/C]ry(-cJe/(C)/C) 

(92) 

As regards the AFMl state in the Coqblin-Schrieffer 
model, its treatment in a general case is a difficult and 
cumbersome problem. However, one can see from (fz^ ) 
that in the leading approximation in 1/N the scaling 
equations coincide with those for the FM state with the 
replacement Jp — s- Jp 



(2) 



8. Improved version of scaling equations: an account 
of dissipative contributions to spin dynamics in or- 
dered phases 

The transition from the dissipative spin dynamics, which 
is characteristic for the PM phase, to the dynamics 
with well-defined spin-wave excitations (ordered FM and 
AFM phases) results in occurrence of singularities in the 
scaling function i]{x) at a; ^ 1 and in a decrease of 
the gc value. One may suppose that in the situation 
of strongly suppressed saturation moment {g is close to 
gc) the character of spin dynamics in the ordered phases 
should change draatically. By the analogy with weak 



itinerant magnetsEj one may expect that for S S a 
considerable part of the localized-spin spectral density 
comes from the branch cut of the spin Green's function 
rather than from the magnon pole. 

In this section we shall demonstrate that this indeed 
takes place provided that our approach is slightly mod- 
ified. To this end we shall analyze the structure of the 
spin Green's function with account of the singular Kondo 
corrections in more detail. 

First we consider the case of a ferromagnet. We obtain 
near the magnon pole 



l/Zq - 1 = - 



\ duj 



(94) 



Besides that, there exists the singular contribution which 
comes from the incoherent (non-pole) part of the spin 
spectral density. To calculate the renormalization of g 
we use, instead of (|40|), the representation of the electron 
self-energy in terms of the spectral density 



J CO q 



E — tk-c 



-Im((5+|5Zq)), 



1 - '^k-q 



Im((5+|5lJ)^ (c 



Thus the magnon pole contribution to g^f is multiplied 
by and the incoherent contribution by 1 — Z. The 
renormalizations of uJef and Z are obtained in a similar 
way. We derive in the nearest-neighbor approximation 



dgcf{C)/dC 
dlnZJef(C)/dC 
d{l/Z)/dC 



-A (96) 
aVA/N (97) 
d\nSef{C)/dC ^VA/N (98) 



where a = 2{\- a), N = 2 

A = [glf{C)/C][Zi^coh{-^ef{C)/C)+{l-Z)i^,^coh{- 

(99) 

'Hcoh = r}^^ , and the function rjincoh is, generally speak- 
ing, unknown. For further estimations we may put 

PM 

'Itncon — '/ 

Although the account of the incoherent part does 
not modify strongly numerical results (see Sect. 10), the 
physical picture of magnetism changes drastically. Ac- 
cording to (98) we have 



ZJef{C)/C)] 



1 



ziO 



1 +ln = 



S 



(100) 



Consequently, the increase of magnetic moment owing to 
the Kondo screening leads to a considerable suppression 
of magnon contributions to the spectral density. Unlike 
the case of weak itinerant magnets, this suppression is 
logarithmic. 

In the case of an antiferromagnet the calculations are 
performed by taking into account the expressions ([5l|), 
(p2|). In the nearest-neighbor approximation (Jq+q = 
— Jq) we obtain 
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— \\coh 



85'^(Jo — Jq, 



(1- 



(101) 



(the quantity (101) just determines tlu 
ization of electron spectrum, cf. Ref.Ej 
for AFM differ from those for FM b 
a — > 1 in the scaling equations 



: Kondo renormal- 
).Then the results 
r the replacement 
I) only. 



^AFSmall dispersion, or damping (the latter occurs in next 
P°Vders in q, cfJ23) results in eliminating this divergence. 
Thus this case needs a more consistent consideration 
with account of higher-order corrections, which will be 
performed elsewhere. 



9. The scaling behavior 
SchriefFer model 



in the large-Af Coqblin- 



It is instructive to consider the limit N — > oo. (To 
avoid misunderstanding, it should be noted that this 
limit with inclusion of spin dynamics in the zeroth ap- 
proximation differs somewhat from the considerations of 
RefsO'E3).Then the renormalizations of spin dynamics 
and Sef are absent, and the transition into the non- 
magnetic Kondo-lattice state cannot be described. How- 
ever, peculiarities of the dependence gef{C) for various 
types of magnetic ordering are described qualitatively, 
explicit analytical expressions being obtained. We have 



l/5e/(C)-l/.9 = G(C)=- 



C' 



(102) 



where D is the cutoff energy defined by gef{—D) = g. 
Performing integration we obtain 



G"^''' {G) 

G^^^\G) 
G^^^'iG) 

G\G)='-\n\{G^-uJ^)lD^\ 



\n\G/D\ 
\n\G/D\ 



-uj/C\ + l 



[G/uj- l)ln|l 

]-[{GVu^-l)\n\l-u^'/G'\ 



e{\G\-uj)M^(\G\ + VC^^)/D) + ^ 



The dependences l/gef{^ — ln|-D/C|) are shown in 
Fig.l. The effective coupling constant gef{C) begins to 
deviate strongly from its one-impurity behavior 



l/g,fiG)^l/g-ln\D/G\ 



(108) 



at |C| ~ uj. One can see that at small uJ <C |C| spin 
dynamics results in a decrease of gef{G) for PM and 
FM cases, and in an increase for AFM case and for the 
case of strong anisotropy. For the dispersionle ss s pin 
dynamics mode the singularity at \G\ = ZJ in (102) is 
non-integrable, and consequently gef{G) diverges at this 
point for arbitrary small g. However, an account of a 



It should be noted that the equation (102) can be used 
even for small N provided that g is considerably smaller 
than the critical val ue q^ (see the next Section) . Besides 
that, the equation (102) works in PM and FM phases 
provided that kp is small, so that, according to (|3^), 
a — > 1. Ho weve r, in the case of the ferromagnet we have 
instead of (|104[) 



G''^^(G)=\n\G/D\- 



N - 1 
N 



(C/cJ- l)ln|l -ZU/CI (109) 



ln((C2 +u2)/^2) + £arctan(^) - 1 



+ — (C/cJ+l)ln|l + cJ/C| + l 

In the 3D AFM case l/gef{G) has a minimum with 
decreasing |C|. For a 2D AFM, 1/gef has a square-root 
singularity at \G\ < ZJ and is constant at |C| > w. As 
we shall see in the next Section, these features retain at 
finite N. The minimum occurs also in the 2D AFM 
case provided that we introduce even small magnetic 
anisotropy. In the 3D case, the anisotropy results in that 
the minimum becomes deeper. For an anisotropic ferro- 
magnet, the minimum is absent at iV — > cxd since only 
the contribution oir]^ survives. However, it is present for 
finite N due to contribution of rj^ . The picture for N = 2 
and g well below gc (when the equation (102) works qual- 
itatively) is shown in Fig. 2. Note that, contrary to the 
case N ^ oo, for FM with N = 2 spin dynamics results 
in a decrease of gef{G) at(il06)too small C. 

The boundary of the strong coupling region (the renor- 
malized Kondo temperatui^ej is determined by G{G = 
-[-jr^)_=g — 1/.9- Of cours^j^Q^ means here only some 
characteristic energy scale extrapolated from high tem- 
^^eri^f^^ff^^^^j^he^d^y^^^^ the ground 

3 a rnore detailed consideration. In the PM 



' state requires 
and FM phases spin dynamics suppresses T^. On the 
other hand, in the AFM case spin dynamics at not large 
lJ results in and increase of T^. 

Provided that the strong coupling regime does not oc- 
cur, i.e. g is smaller than the critical value gc, gef{G 
0) tends to a finite value g*. To leading order in \n{D/uj) 
we have 



l/.g, = A = lii{D/U) 



(110) 



However, an account of next-order terms results in an 
appreciable dependence on the type of magnetic ordering 
and space dimensionality. For PM, FM and 2D AFM 
phases the critical value gc is given by 1/gc — ~G(0). 
Then we obtain 
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l/g* = l/g-l/5c 



with 



l/5c = A + 



1 PM, FM 
In 2 2D AFM 



(111) 



(112) 



At the same time, in the 3D AFM case (and also in the 
anisotropic 2D case) gc is determined by the minimum of 
the fmiction G{C). Thus g* remains finite dA, g gc — 0- 



l/g* = l/g - l/g, + l/g* 

where 

l/.gc = A+(l + ln2)/2, 
l/g,* = A = (ln2)/2~0.35 

A being the depth of the minimum. 
For C ^ we have 



(113) 
(114) 



G(C) - G(0) 



'7r|C|/2cJ PM 
\C/uJ\\n\U/C\ FM 
-{C/ujf In \uj/C\ 3D AFM 



(115) 



(however, G^^(C) - G^*^(0) ~ (C/ZJ)^ for the function 



( |109D at TV = 2). As follows from ( |115| ), for 5 ^ 5, + 
we have in PM and FM cases to logarithmic accuracy 



-9 -9c 



K 



OJ- 



9t 



2/7r, PM 
H{g~g,)/g,], FM 



(116) 



On the other hand, in the AFM phases is finite at 
g gc + 0: 



K 



LO X 



1, 2D AFM 
l/%/2, 3D AFM 



(117) 



10. Scaling behavior for finite A*' 

Some results of the large-A^ limit hold in the finite-A^ 
case too, but there occur a number of new important 
features. To consider the general case we write down 
the first integral of the system (89), (|90| ) 



gef{C) + {N/Va)\ni:j,f{C) = const 
which results in 

W,f{C) = LoeM-{aV/N)[g,f{C) - g]) 
As follows from (|^), 1^ 



(118) 



(119) 



Sef{C) 



ef{C) 



1/a 



(120) 



Substituting ( |119| ) into (|89|) we obtain 

9(l/5e/)/9^ = -^{X + iaV/N)[g,f ~ g] ~ (121) 
where 



*(^) = f]{e-^), e = In \D/C\, A = \n{D/uJ) > 1 

Prese nce of g,/ in the argument of the function ^' in 
(121) leads to drastic changes in the scaling behavior in 
comparison with the large- limit. We describe below 
the renormalization process in various cases. 

If g is sufficiently large, the increase of gef will lead 
to that the argument of 4" will never be small, so that 
the scaling behavior will be essentially the one-impurity 
one. At smaller g the function l/gef{£,) begins to devi- 
ate strongly from the one- impurity behavior l/gef{C) = 
l/g — £, starting from ^ ~ ^1 where ^1 is the minimal 
solution to the equation 



X+{aV/N)[g,f{0~g]^C 



(122) 



If the argument remains negative with further increasing 
^, 1/gefiO will decrease tending to the finite value g*, 
the derivative d{l/ gef) /d^ being exponentially small, so 
that the situation is close to the large- A^ case (Fig.l). 
However, {aV/N)g can increase more rapidly than 
^, so that the second solution to (122), ^2, will occur, and 
the argument of the function 5* becomes positive again. 
Then 5e/(C) will diverge at some point ^* = ln|£>/r^|. 
The divergence is described, unlike the large- A^ limit, by 
the law 



ef{0^i/{s,*-i) 



(123) 



since ri{x <C 1) = 1. The behavior ( |123| ) takes place 
starting from ^ ~ ^2- 

The dependences gefiO the PM case at small ji?— (/d 
are shown in Fig. 3. The behavior gefiS,) between ^1 and 
^2 is nearly linear, but is somewhat smeared since VE'(^) 
differs considerably from the asymptotic values and 1 
in a rather large interval of ^. 

The case of magnetically ordered phases has a number 
of peculiarities. Here the singularity of the function 5'(^) 
at ^ = turns out to play the crucial role. In particular, 
one can prove that g^f diverges at some ^ at arbitrar- 
ily small g (i.e. g, ~ 0) unless the singularity cutoff is 
introduced. Indeed, when approaching the singularity 
point with increasing ^, the derivative dgef/d^ rapidly 
increases, and the argument of the function ^> in (121) 
inevitably starts to increase at some point ^1. Thus the 
singularity point cannot be crossed, and the argument of 
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the function 5* is always positive. Since ^'(^' > 0) > 1, 
we have d{l/gef)/d^ < —1 at arbitrary Therefore the 
effective coupUng constant diverges at ^ = ^* < 1/g. 

To make the value of Qc finite one has to cut the sin- 
gularity of the scaling functions. This may be performed 
by introducing small imaginary parts, i.e. replacing in 
l2|) 



{aV/N){dgef/da^,^ = (aW^^).gi*max = 1 (127) 

On the other hand, we may estimate from the linear 
asymptotics l/gef{£,) — I/5 — ^ (which holds up to ~ 
(,~\+iaV/N)g,) 



1/g, ~ 1/g - a 



(128) 



lull 



(1 



^'4 

-1/20(1 



ln[(l 



Comparing (|128| ) and (127) we derive the rough estima- 
(124) tion 



Re(l 



iS) 



-1/2 



= {[(1 - xf + <52]i/2 + 1 _ a;]/2}i/2/[(l - + S^]'/^ 



Then ^ becomes bounded from above: 



^(1) 



iln^ FM 



In (5 



3D AFM 



i(5-i/2 2D AFM 



(125) 



Other way to cut the divergence in the 2D AFM case is 
to retain strictly the proportionality to the step function 



(l-x) 



-1/2/ 



[l-x)^Re{l-x + iS) 



-1/2/ 



[l-x) 
(126) 



so that ^(^) cannot take small values. 

The value of 6 should be determined by the magnon 
damping at g |k-k'| ~ 2kF (see (Q), (|l|)). This 
damping is due to both exchange and relativistic inter- 
actions. The damping owing to exchange scattering by 
conduction electrons should be formally neglected within 
the one-loop approximation, since this contains a more 
high power of /. The magnon-magnon interaction in the 
Heisenberg model gives the damping at non-zero tem- 
peratures only. However, the relativistic (e.g., dipole) 
interactions give a damping at T = due to zero-point 
oscillations. Hereafter we put in numerical calculations 
5 = 1/100. 



The behavior of the solutions to (121) in magnetic 
phases for g well below gc is similar to that in the 
large-iV limit (see Fig. 2). In particular, the function 
l/gef{£) has a minimum in the 3D AFM case with 
the same depth, and a similar situation occurs in the 
presence of anisotropy. The presence of the minimum 
may result in non-monotonic temperature dependences 
of physical quantities which are sensitive to the Kondo 
screening, e.g., of the effective magnetic moment. These 
dependences are o btain e d qu alitatively by the replace- 
ment |C| T in ([TT9|),([l20|). Of course, the standard 
monotonic spin- wave corrections should be added to the 
Kondo contributions. 

It is important that, except for a very narrow region 
near g^^ at approaching g^ the value of g* becomes prac- 
tically constant, g* ~ g(Ci) = 5i- This value is estimated 
as 



1/g, = A + {aV/N)g, + 1/g, 
= A+(iV*„,ax/al^)"'/' 



(129) 



The value of 1/ gi can be also estimated from the linear 
asymptotics (123): 



17.92 ~ r - 6 



(130) 



Owing to the singularity, ^2 is close to ^1 except for very 
small \g — gd- In the latter case the behavior g(^i < ^ < 
^2) is practically linear. 



(ay/Ar)<?e/(a <C<6)-e-A 



(131) 



(this behavior is discussed in detail in the next sections). 
This is explained by that the argument of the function 
'J in (121) shoul d be nearly zero. The behavior at ^ > ^2 
is described by (123). We may estimate from (130) 



r ~ 1/g + l/<72 - 1/51 



(132) 



For the cutoff ( p^ ) with (a; > 1) = in the 2D AFM 
case (and in similar situations) ^* tends to a finite limit 
^* at (7 ^ (7c + 0. Indeed, as follows from the form of 
the scaling function (82), one has up to the divergence 
point 4'(^) > 1, so that i* < 1/g < 1/gc- In such a 
situation g* turns out to be also finite at 5 — > ~ 0, 
and the character of approaching g, is quite different 
from that in the large-A^ limit. With increasing the 
position of the l/gef{S/) singularity point is shifted to 
right due to ra pid in crease of g^f in the argument of the 
function * in ( |l2l|) . The shift should stop at ^1 < C- 
This takes place just at g — > .9c — 0. Thus l/gef{C} should 
vanish discontinuously at g = g,. The "maximum" value 
of ge/(Oi 9eT^ = 9c^ is estimated from 



A+(ay/iV)(5;7''-5c)-Ci=0 



(133) 



Since for g > g, th e de crease l/(7e/(0 C > Ci is prac- 
tically linear (see (123)), we may estimate 



C - Ci ^ 1/5, 



max 



(134) 



If we accept the cutoff (124), a small increase of ^* and 
1/g* will take place in an extremely narrow region near 
gc- This increase cannot be practically observed. 
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In the 3D AFM case the statement about the finiteness 
of ^* does not, strictly speaking, hold since there exists 
a small region where ^ (^) is positive and takes arbitrar- 
ily small positive values, so that the decrease of l/gef{0 
can be slow. However, the behavior £,*{g) is in fact de- 
termined by the logarithmic singularity of the function 
^'(^) except for a very narrow region near gc- The nu- 
merical calculations yield the estimation g — gc ^ 10^'* 
for the region where ^* starts to increase. Thus, from 
the practical point of view, we may put ^* to be finite. 
The corresponding value of g"^^^ = gi is determined by 
( |133D and gl is smaller than g^f'^. 

At very small g — gc ^ 10~^ the value of g* starts to 
increase. However, due to the minimum of the function 
l/.9e/(C)j 9c ~ remains finite at 5 — > gc, as well as 
in the limit iV — > cxd. A similar situation takes place for 
anisotropic ferromagnets. 

In the case of an isotropic ferromagnet the influence 
of the singularity is somewhat weaker since ^'(^) does 
not change its sign and takes small positive values up 
to infinity, so that ^* starts to increase appreciably at 
5-5,-10-3. 

The dependences gefiO magnetic phases according 
to (121) at g ^ gc i are shown in Figs. 4-6. 

Now we consider the results of the approach of Sect. 8, 
which takes into account the incoherent part of the spin 
spectral density The integral of motion of th e system 
of equations (pq), ( |9^ ) has the same form (118), and we 
obtain 



{V/N){^ - X) ^ZT 



(138) 



1 _ 1 LU 



l + ^[5e/(e)-5] (135) 



The corresponding equation for gcf reads 
dil/9ef) 



- *»„co/i(A + {aV/N)[gcf ~ g] - 



-[^cok{X+(aV/N)[gcf-g]-0 



V 



As demonstrate numerical calculations, the account of 
incoherent contribution results in a smearing of the non- 
monotonous behavior of gefiO in the 3D AFM and 
anisotropic cases, and in some region of I5 — 5c| the mini- 
mum of 1/ gcf (0 vanishes completely. Therefore — > 
at g gc-, unlike the situation for the equation (121). 

The influence of the incoherent contribution on ^* 
and g* is considerably suppressed by the singularity of 
the function "^cohiS,)- The region where this contribution 
starts to play a role is determined by the quantity S. In 
particular, for the 2D AFM case its influence on ^* is 
practically absent since the divergence of gefiO occurs 
due to the singularity of 'i'cohiO- Note that since g* is 
finite at g gc, the coherent contribution survives up 
to 5c- 

The comparison of the results of various approxima- 
tions is presented in the Table 1. One can see from this 
Table that for = 2 the relation of the gc values in the 
ordered phases an in the PM case is reversed in com- 
parison with the limit N 00. This fact is due to the 
influence of the scaling function singularities. It should 
be noted that at larger 5 ~ 1/5 the value of gc in the 
ordered phases exceeds 5^*^, as well as in the large- TV 
case. In the case (b) the value of gc is intermediate be- 
tween gc"'' and 5^*^ and closer to gc'^\ With increasing 

a or A^ the difference between gi''^ and gi'^^ becomes still 
smaller. 

Table 1. The critical values gc and ^* for d i fferen t 
i nagn etic phases in the cases N = 00 (see ( |112| ), (114), 



(117)) and A^ = 2 in the approximation of Sect. 7 (a) and 
with account of the incoherent contribution (b). The 
parameter values arc A = 5, a = 1/2, S = 1/100. For 
(1/Sr6)= 2, the "critical value" of Q is estimated from the 
plateau in the dependence S,*{g) (see the discussion in 
the text). 



-^^ncohiX + iaV/N)[gcf - 5] - 0]/[l + ^(5e/ 

The role of the incoherent contribution becomes impor- 
tant only provided that Z deviates appreciably from 
unity, i.e. gcf — g is large. This takes place in a rather 
narrow region of I5 — 5c|- The estimation for gi in the 
case of small 6 reads now 







PM 


FM 


3D AFM 


2D AFM 


^Jv ^ 00 


9c 


0.167 


0.167 


0.171 


0.176 


Cc 






5.35 


5 


A^ = 2 (a) 


9c 


0.154 


0.139 


0.132 


0.127 


Cc 




6.13 


6.07 


6.07 


A^ = 2 (b) 


9c 


0.154 


0.141 


0.136 


0.131 


Cc 




6.23 


6.17 


6.16 



aV9l/N 
1 + Vgi/N 



^ con ^ 



(137) 



so that, as follows from ( |129| ), gc increases. The de- 
pendences gefiO 3't small \g — gd according to ( |136| ) 
are shown in Fig. 7. One can see that a crossover from 
the well-linear "magnetic" behavior to a PM-like "quasi- 
linear" behavior occurs with increasing ^. The point of 
the crossover is estimated from 2' 4'™^'^ ~ 1, i.e. 



The dependences t/g*{g) and £,*{g) according to ( |136| ) 
are shown in Figs. 8- 11 (of course, these Figures do not 
show the above-discussed increase of \/g* and ^* in the 
AFM case, which takes place at very small I5 — 5c|)- The 
experimentally observable quantities can be obtained 
from these data by using the formulas 

S*=Scf{C = 0) = ScM-iV/N)[g* -g]) (139) 
W* = ZJcfiC - 0) = LJcxp{-(aV/N)[g* - g]) 
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id < 9c)- For g > gc we have 



so that 



(140) 



-fiN/Va)Hg,~g) 



(144) 



One can see from Figs. 8-11 that, provided that g is far 
from gc, we have the one-impurity behavi or C* (g ) — Vffi 
and the dependence l/g*{g) is given by ( |111| ), ( |ll3| ), as 
weh as in the large- iV hmit. 



11. Critical behavior on the boundary of the strong- 
coupling regime: breakdown of the Fermi-liquid pic- 
ture 

To investigate the "critical" behavior of ^* at g — > gc we 
consider the function l/gef{g,C) at small Ig — gd, |C|. 
First we consider the results of the solution of the equa- 
tion (121) . W ith approaching gc, the ^-region, where the 
behavior (123) takes place, becomes very narrow and not 
too important for determining Unlike the exponen- 
tial (in ^) behavior in the large-iV limit, we have from 



gef{g^gc,^^^)^-{N/Va)^ 



(141) 



However, the dependence ^*{g) turn s ou t to be qualita- 
tively the same as for ^ oo (see (116)), 



r - 7 - 9c) 



(142) 



Numerical calculations yield 7=1/2 for FM at TV = 2, 
and 7 = 1 for FM with > 2 and FM. These values 
are the sa me a s for a ^ (or according to the large- 
equat ion ( |102[ ), provided that we take for FM the func- 
tion ( |109| )). Thus one may put forward the hypothesis 
that the critical exponents are universal, i.e. depend on 
the type of magnetic ordering only, but not on N, V and 
a. In the AF M p hases numerical calculations yield the 
dependences (142) with 7 ~ 0.1 {d ~ 3) and 7 ~ 10~^ 
(d = 2). 

At the same time, the behavior of g* sd, g —i- gc — 
changes in comparison with the large- iV limit. It turns 
out that for finite N one may establish a scaling relation 
of relevant variables at g > gc and g < gc, &s well as in 
the standard theory of critical phenomena. To find this 
relation, we consider our problem in the region Ig — gd > 
e where e — > determines a scale of approaching to 
the critical point. When crossing the cut region, the 
argument of the function 5* should not shift considerably 
(Figs. 3-6). Indeed, this argument must be close to zero; 
in the ordered phases it is fixed by the singularity point, 
and for PM a considerable "smearing" takes place. Then 
we may estimate 



As demonstrate numerical calculations (Fig. 8), for the 
PM phase the increase of I/5* at 5 not too close to gc is 
almost lin ear in g, as well as in the large A^- limit, and the 
behavior (144) takes place only starting from g* ^ 10 
which is, strictly speaking, beyond the applicability of 
one-loop scaling. At the same time, for the FM case 
the logarithmic dependence takes place in a considerable 
interval oil/g* . As discussed in previous Section, for 2D 
AFM the increase of g* and ^* at g ^ gc can be hardly 
observed. 

Of course, f or 3 D AFM 5* = 1/A is in fact finite, and 
the behavior ( 144 ) takes place at not too large g* . More 
exactly, we can write down 



l/g* c^l/YtiN/Va) ln(.g, - g)] + A 



(145) 



However, practically the "saturation" region is extremely 
narrow and cannot be achieved because of the smallness 
of 7. 

When taking into account the incoherent contribution, 
a crossover to a PM-like regime takes place aX g ^ gc in 
the dependences g*{g) and ^* (g), so that at very small 
I5 — 9c\ we have the beh avior (]142|), (144) with 7 = 1. 



Basing on the results (142) 
sume that at g ^ gc, C 
scaling behavior 



144), it is natural to as- 
— one has the one-loop 



gcf{g, C) = -(A^7/aV') ln(|C/cJ|i/^ + B{gc - g)/g) 

(146) 

(B > is a constant, the argument of the logarithm 
should be positive). Then, according to (119), (020) 



Zocfig, C)lu5 ^ (|C/cJ|i/^ + B{gc - g)/g)-' (147) 

Scf{g,C)/S ^ {\C/ZJ\'^^ + B{gc~ g)/g)-^/'' (148) 

In particular, we have at g — s- g^ the power-law depen- 
dences 



T^,uj* - [±{g-gc)]"' 
S* - {9c - 9)'''" 



(149) 
(150) 



Thus the "critical exponents" for the characteristic en- 
ergy scales, na mely, T'^ at g > gc and TD* at g < gc 
coincide. Using ( |l48| ) we obtain at g = g^ C ^ Q 



Scf{C) _ fZJcfiC) 



X+{aV/N)[g*-g]:^C 



(143) 



1/a 



\c\ 



lla. 



(151) 
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One of the most interesting consequences of our pic- 
ture is a possibility of a non-Fermi-liquid behavior on 
the boundary of the strong coupling region. Indeed, dur- 
ing the renormalization process at C the effective 
spin-fluctuation frequency tends to zero, and the corre- 
sponding spectral density is concentrated near ui — 0. 
In this sense, the situation is close to that in the one- 
impurity two-channel Kondo problem where a collective 
mode with zero frequency occurs, which leads to a break- 
down of the Fermi-liquid picture due to electron scatter- 
ing by this "ultrasoft" modc£3. Unfortunately, our per- 
turbation approach does not permit to determine explic- 
itly the temperature dependences of observables since 
the coupling constant is not small in this regime. How- 
ever, the calculations may be performed within the large- 
[l] s — f model (see the next Section). 

Of course, vanishing of uJef{C) at |C| ~ Tj*. and of 
and W* at g = Qc is the result of the one-loop scal- 
ing, i.e. of using the lowest-order perturbation theory 
at derivation of the renormalization group equations. 
In fact, one may e.Yuect that in the strong-coupling re- 
gion ZUe/(C) ~ TkEj'U. One may assume that the cor- 
rect scaling behavior, which may be continued into the 
strong-coupling region, differs from the one-loop behav- 
ior by the replacement 



(152) 



A scaling law, which is more general than (147), could 
be expected to have the form 



cJe/(c)/cu = (r^/cj)0(c/r^) 

Then one has 

S,f{C)/S = 4,{Lj,f{C)lLj) 



(153) 



(154) 



A detailed investigation of magnetic properties (in par- 
ticular, of the formation of small moments, which are 
characteristic for heavy-fermion systems) reduces to de- 
termining an explicit form of the functions <j) and 
This problem cannot be solved within perturbative ap- 
proaches. 



12. The non-Fermi-liquid behavior in the degenerate 
s — / model 

As we have seen in Sect. 9, in the large- iV limit the renor- 
malization of magnetic characteristics is weak in compar- 
ison with that of electron spectrum. An opposite situ- 
ation occurs in the case of large [I] where the number 
of electron branches is much larger than that of spin- 
wave modes, so that the renormalization of spin dynam- 
ics plays the crucial role. 



It is instructive to consider the large-Z limit in the s — / 
model {N = 2) with 



[I] 



D, <? ^ 0, [%V2 = ff" = const 



Then the effective s — f interaction is unrenormalized, 
Ijef = = const, and the scaling equation takes the form 



dx 
where 



ag'^iX + X-0 



(155) 



x(e) = In 



When taking into account the incoherent contribution 
we obtain instead of ( |155| ) 

^ = ag^[Z-^,oh{\ + X - + (1 - ^)«'mco/.(A + X - 0] (156) 
Z=l/{l + x/a) 
By introducing the function 

UJef 



(157) 



the equation (155) can be readily integrated to obtain 



, l-a52*(e) 



(158) 



However, a simple qualitative analysis can be performed 
immediately for both the equations (155) and (156). 
In the PM phase we have 



X(0 ^ ag'i 
up to the point 

1 - ag^ 

Thus a power-law behavior occurs 
cJe/(C)c.cJ(|C|/Z?)^ /3 = a5?2 

For i > a, 

x(0^x(ei) = Aa5^(l-ag2)-i 



(159) 



(160) 



(161) 



(162) 



is practically constant. 

Note that unlike the case I — 0, which was discussed 
in the previous Section, the bare coupling constant g and 
the exponent /3 can be sufflciently large. At ^ ~ 1, ZJef (C) 
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decreases rather rapidly during the renormahzation pro- 
cess. Thus a "soft mode" situation occurs, which may 
lead to a NFL behavior. 

To investigate modification of the electron spectrum 
we may calculate the second-order perturbation theory 
corrections, which are formally small in Replacing 
uJ uJef{C) in the usuaLsecond-order result for the elec- 
tron self-energy (cf. RefBj) and introducing the effective 
electron density of states at the Fermi level Nef{C) we 
obtain 



NefiC) 



ZjC) 

ZJefiC) 



D_\' D_ 
\C\) 



(163) 



To investigate qualitatively temperature dependences we 
may replace |C| — > T. Thus one may expect an essentially 
NFL behavior of the electronic specific heat, 

Ce{T) ^ TNef{\C\ ^ T) ~ T^-f'h-iiD/T), (164) 



magnetic susceptibility, 

XUT) - Slf{\C\ -> T)/T ^ \IT^~^I\ 



(165) 



transport properties etc. 

In magnetically ordered phases, the situation for ^ > 
^1 changes since the singularity of ^'co/i(0 pla^ys an im- 
portant role. Provided that dip'^^^^ > 1, the argument 
of the function ^coft, at ^ > becomes almost constant, 
and we obtain 



(166) 



The behavior ( 166 ) is similar to the dependence (|147|) for 
finite I, and corresponds to g = gc- In the case of equation 
(155), this behavior takes place up to ^ = oo. On the 
other hand, an account of the incoherent contribution 
results in that the increase of x stops at ag^'i''^^^ — 
1/Z = 1 + x/cL, i.e. at 



" coh 



?2 = A + Xr 



(167) 



Thus the value of ^2 is determined by the quantity 5. 
The dependence x(C) for a 2D antiferromagnet is shown 
in Fig. 12. In the presence of the in coher ent contribution 
the region, where the dependence (166) holds, is rather 
narrow (especially for not too small S). However, a more 
exact consideration of spin dynamics (rather than using 
the spin diffusion approximation) may change consider- 
ably the results. 

In the regime (166) we have the result ( |163D with P — 1 
and D ^uj, so that in the AFM case (a — 1) one obtains 



Xm{T) ^ const, Ce(T) - ln(cJ/T) 



(168) 



(169) 



However, the NFL behavior takes place in a wide tem- 
perature region Tk <C T < ZZJ. At T < Tk the renor- 
mahzation of gef becomes important and, as discussed 
in the end of the previous Section, a more complicated 
scaling behavior may take place. 

As discussed in Section 11, a NFL behavior takes place 
even ior I — for g ~ gc- The NFL region becomes 
broader with increasing [I]. The dependences <7e/(C) for 
I = 3 are shown in Fig. 13. One can see that the linear 
dependence with a small coefficient takes place up to ^ ~ 
5, then this is changed by the linear "coherent" behavior 
which is further smeared by the incoherent contribution. 



13. Conclusions 

In conclusion, we resume main results of our consider- 
ation and their relation to properties of the anomalous 
/-systems, and discuss some unsolved problems. 

Three regimes are possible at / < depending on the 
relation between the one-impurity Kondo temperature 
Tfc and the bare spin-fluctuation frequency ZJ: 

(i) the strong coupling regime with Ief{C ^ 0) 00 
where all the conduction electrons are bound into singlet 
states and spin dynamics is suppressed. This regime is 
expected to occur provided that cJ ^ Tk- 

(ii) the regime of a "Kondo" magnet with an apprecia- 
ble, but not total compensation of magnetic moments, 
which corresponds to small \g — gc\- 

(iii) the regime of "usual" magnets with small loga- 
rithmic corrections to the ground state moment andilJe/ . 
(Note that the same situation takes place at / > Ollj.) 

The formation of magnetic state takes place at 



K 



D exp(- 



l/5e) ^^^[^ 



l/Oc-l 



(170) 



(flc = Xgc, A is of order of unity). For N ^ 00 we have 
flc = 1-|-0(1/A) and the strong coupling region boundary 
is determined by the condition Tk — AZJ. For finite N w e 
always have Oc < 1 (see Table 1) and, according to (170), 



For finite, but large [I] the picture discussed fails below 



T^ <^ UJ. Numerical solution of the scaling equations for 
the case N = 2 demonstrates a considerable dependence 
of the critical value gc on the type of magnetic ordering, 
space dimensionality and the structure of the magnon 
spectrum (presence of the gap). 

11 should be stressed that the Doniach criterion gc — 
0.40, which was obtained for a very special case in a sim- 
plified one-dimensional model, cannot in fact be used for 
real systems since gc turn out to be sensitive to param- 
eters of exchange interactions, type of magnetic order- 
ing, space dimensionality, degeneracy factors, magnetic 
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anisotropy etc. It is worthwhile to note in this connec- 
tion that it is a common practice to treat the interplay 
of spin dynamics and the Kondo_effect within the two- 
impurity problem (see, e.g., Refs.EZi). At the same time, 
we have demonstrated that the most important features 
of the scaling behavior are connected with peculiarities of 
the spin spectral density and are not described by this 
model (where the spectral density is a delta-like peak 
corresponding to singlet-triplet transitions). 

We have used in our calculations the simplest "De- 
bye" approximation for the magnon spectrum, i.e. the 
long-wave dispersion law in the whole Brillouin zone. 
Competing interspin interactions owing to the oscillat- 
ing behavior of the RKKY exchange, which may lead 
to frustrations, might be also important for explaining 
magnetic structures in "usual" Kondo lattices. 

The effective Kondo temperature determines a 
characteristic energy scale of the "heavy-fermion" be- 
havior at low temperatures. This may differ consider- 
ably from the one-impurity value Tk, so that in the PM 
phase < Tk, and in the magnetically ordered phases, 
except for very small g — gc, > Tk (see Figs. 9-11). 
In the AFM case depends weakly on g and therefore 
on Tk, and is determined by in a wide interval of g. 

Although our consideration was performed for T = 0, 
one may expect by the analogy with the one-impurity 
problem that the dependences S{T) (in the PM phase, 
S{T) is the local moment determined from the magnetic 
susceptibility) and uj{T) may be qualitatively obtained 
by the simple replacement |C| —> T. 

For N = 2 and small enough 6 (spin excitation damp- 
ing) the Kondo screening in the AFM and FM phases 
is stronger than in the PM phase. This results in an 
increase of the Rhodes- Wohlfarth ratio (ratio of the ef- 
fective moment, as determined from the Curie constant, 
to the saturation moment). For example, if the bare cou- 
pling constant is close to its critical value in the ordered 
phase and correspondingly lower than in the PM phase, 
the ground state moment may be small in comparison 
with the high-temperature one (a behavior, typical for 
most Kondo magnets, as well as weak itinerant electron 
d-magnets). The same situation takes place for the char- 
acteristic spin-fluctuation frequency ujef- It should be 
noted that an increase of spin-fluctuation energies with 
temperature is indeed obseriied in a number of anoma- 
lous /-systems, e.g., U2Zni7E£l. 

Presence of the factor 1 /N in the AFM2 case or large 
value of 6 (magnon damping) may result in that the sup- 
pression of the effective moment and spin-dynamics fre- 
quency turns out to be weaker than in the PM phase. 
It looks like the decrease of moment at magnetic dis- 
ordering with increasing temperature. Such a decrease 
is typical for strong itinerant magnets where it is due 
to the change of electron spectrum at disordering (see 



RefsH3'E3). In the case under consideration this phe- 
nomenon has a quite different (essentially many-particle) 
nature. 

Near the boundary of the strong-coupling region {g 
gc) the relevant variables demonstrate a non-trivial scal- 
ing behavior as functions of I5 — 5c| and ^. In particular, 



TKig ^ 9c + 0), w*(.g gc - 0) - |g - 



(171) 



the exponent 7 depending on the type of magnetic or- 
dering. These results may be of interest for the general 
theory of metallic magnetism. The description of the 
state with small magnetic moments {g gc) turn out 
to differ considerabUi-from that in the theory of weak 
itinerant magnetismE^I. It is interesting that the "crit- 
ical exponents" in the dependences of the moment on 
the coupling constant ( |150D and C (151) turn out to 
be non-integer. The corresponding dependences S{T) = 
S'e/(|C| T) describe an analogue of the "temperature- 
induced magnetism"E3. 

As mentioned in the Introduction, high sensitivity 
of the magnetic state in heavy-fermion systems to ex- 
ternal factors is explained by that in the case (ii) the 
magnetic moment changes strongly at small variations 
of the bare coupling constant. According to our con- 
sideration, the regime of small magnetic moments oc- 
curs in a very narrow region of bare parameters only. 
Provided that g is not too close to gc, a characteris- 
tic interval of the change of the quantity g* (which de- 
termines the renormalized values of magnetic moment 
and spin-fluctuation frequency) by unity is estimated as 
^9^9^^ 9- In the immediate vicinity of gc (where 
the behavior g* ^ —\Ti\g — gc\ takes place) this interval 
becomes still more narrow; 5g ^ \g — gc\. A more con- 
sistent treatment of this regime with account of possible 
renormalization of the scale I5 — 5c| itself requires using 
more complicated (e.g., numerical) scaling approaches. 

"Softening" of the spin excitation spectrum in the crit- 
ical region may result in a non-Fermi-liquid behavior. 
Although the NFL state itself cannot be described within 
the framework of our perturbative approach, the con- 
clusion about the NFL behavior near the boundary of 
magnetic [g < gc) and non-magnetic {g > gc) phases 
seems to be important. This conclusion is confirmed 
by the fact that a violation of the Fermi-liquid picture 
in anomalous /-system is—ceally observed near the on- 
set of magnetic orderingl3. It is difficult to explain 
this fact within the fccquently used one-impurity two- 
channel Kondo modeO. At the same time, our scenario 
of the NFL state formation takes into account essentially 
many — center nature of the system. The width of the 
region where the NFL behavior occurs increases with 
increasing the degeneracy factor [I] and decreases with 
increasing A''. In the formal limit [I] — > 00 the pertur- 
bation theory in g remains applicable, so that explicit 
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expressions for thermodynamic and magnetic properties 
can be obtained. 

On the whole, the physical behavior, which occurs as 
a result of an interplay between the Kondo effect and 
intersite exchange interactions, turns out to be very rich 
and differs for various model versions. 

The work was supported in part by Grant 96-02-16000 
from the Russian Basic Research Foundation. 

Figure captions 

Fig.l. The dependence l/(?e/ on ^ = ln|D/C| in the 
large-iV limit with A = ln(£'/ZZ;) — 5, g — 0.15 for a 
paramagnet (dashed line) and different magnetic phases 
(solid lines): (1) ferromagnet (2) 3D antiferromagnet (3) 
2D antiferromagnet (4) magnet with a strong anisotropy 
LUq — > oo. 

Fig. 2. The dependence l/gef{£,) according to the 
equation (121) at = 2 in the anisotropic case (wq = 
0.2) for different magnetic phases: (1) ferromagnet (2) 
3D antiferromagnet (3) 2D antiferromagnet, and for the 
isotropic ferromagnet (dotted line) and 3D antiferro- 
magnet (dashed line). The bare coupling parameter is 
g = 0.11, other parameters are the same as in Sect. 10. 

Fig. 3. The scalin g trajectories gefiO in a paramagnet 
according to (^2l]) with iV = 2 for .g = 0.153917 > g^ 
(upper line) and g — 0.153916 < gc (lower line). 

Fig. 4. The s calin g trajectories gef{£,) in a ferromagnet 
according to ( |l2l| ) with iV = 2 for g = 0.13868 > gc 
(upper line) and g ~ 0.13867 < gc (fower line). 

Fig. 5. The scaling traje ctori es gefiO in a 3D anti- 
ferromagnet according to (121) with iV = 2 for g = 
0.1320382 > gc and g = 0.1320381 < gc . 

Fig. 6. The scaling traje ctori es gefiO in a 2D anti- 
ferromagnet according to (121) with N — 2 ioi- g — 
0.1266714 > gc and g = 0.1266713 < gc- 

Fig. 7. The scaling traje ctori es gefiO in a 3D anti- 
ferromagnet according to (136) with N = 2 ioi g = 
0.136305 > gc (upper line) and g — 0.136305 < gc (lower 
line). 

Fig. 8. The dependences 1/ g*{g) for g < gc and £,*{g) — 
X for g > gc in a, paramagnet (A = 5). The dashed line 
is the curve 1/g — A. 

Fig. 9. The dependences l/g*{g) for g < gc a nd ^* (g) — 
X for 5 > in a ferromagnet according to (136) The 
parameters are X = 5, S = 1/100. 

Fig.lO. The dependences l/g*{g) for g < gc and 
^*{g) — X f or g > in a 3D antiferromagnet accord- 
ing to 

Fig.ll. The dependences l/g*{g) for g < gc and 
^*{g) — X f or g > (7c in a 2D antiferromagnet accord- 
ing to (III). 

Fig. 12. The dependence of x = ln(w/we/) vs. ^ = 
ln|Z?/C| for a 2D antiferromagnet at [I] = oo with S — 
lO^'^, X — 5, g — 0.6) according to (|155|) (dashed line) 



and with account of the incoherent contribution (solid 
hue). 

Fig. 13. The scaling trajectories gefiO in a 3D an- 
tiferromagnet according to (136) with N — 2, I = 3, 
S ~ 1/100 for g = 0.10556 > gc (upper line) and 
g — 0.10555 < gc (lower hue). 
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